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Abstract. We study the mapping from Lagrangian to Eulerian space in the context of the Effective 
Field Theory (EFT) of Large Scale Structure. We compute Lagrangian displacements with Lagrangian 
Perturbation Theory (LPT) and perform the full non-perturbative transformation from displacement 
to density. When expanded up to a given order, this transformation reproduces the standard Eulerian 
Perturbation Theory (SPT) at the same order. However, the full transformation from displacement 
to density also includes higher order terms. These terms explicitly resum long wavelength motions, 
thus making the resulting density field better correlated with the true non-linear density field. As a 
result, the regime of validity of this approach is expected to extend that of the Eulerian EFT, and 
match that of the IR-resummed Eulerian EFT. This approach thus effectively enables a test of the IR- 
resummed EFT at the field level. We estimate the size of stochastic, non-perturbative contributions 
to the matter density power spectrum. We find that in our highest order calculation, at redshift 
z = 0 the power spectrum of the density field is reproduced with an accuracy of 1% (10%) up to 
k = 0.25 ft,Mpc~^ {k = 0.46 fiMpc”^). We believe that the dominant source of the remaining error is 
the stochastic contribution. Unfortunately, on these scales the stochastic term does not yet scale as 
k* as it does in the very low k regime. Thus, modeling this contribution might be challenging. 
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1 Introduction 

The Large Scale Structure, i.e., the distribution of matter and galaxies on large scales, has the 
potential to constrain the history and composition of the Universe in a way that is complementary to 
the Cosmic Microwave Background (CMB). It even has the potential to tighten CMB constraints on 
the physics of the early and late accelerated expansion. The latter can be explored through precise 
measurements of the expansion history using the Baryon Acoustic Oscillation (BAO) [1] method, while 
inflation and the generation of the seeds for structure formation can be constrained using primordial 
non-Gaussianities [2, 3]. In order to be able to extract the full potential of present and upcoming 
surveys in answering these fundamental questions, we require accurate predictions for the clustering 
statistics. On the largest scales linear theory correctly describes the growth of structure, but on 
smaller scales, non-linear effects become important before on the smallest scales overdensities collapse 
to form virialized objects. Numerical simulations of structure formation have become a powerful 
tool to model non-linearities and interpolation techniques have been developed to allow for a fast 
exploration of the parameter space [4, 5]. Yet, numerical simulations hide some of the physics and 
convergence between codes and parameter choices struggle to pass the percent barrier [6, 7]. 

Perturbative techniques (for a review see [8]) are an alternative, extending the validity of linear 
theory into the weakly non-linear regime and improving precision on large scales. Recent years have 
seen a resurgence of interest in perturbative approaches to study the development of structure in our 
Universe. A new development has been the introduction of Effective Theory techniques [9-14] and a 
substantial effort went into higher order computations and testing them against numerical simulations 
[15-18]. These approaches have the potential of providing very accurate calculations for observables 
on large, linear and mildly non-linear scales. This work follows our recent paper [18] (hereafter BSZ), 
where we studied the perturbative solution for the displacement held of dark matter particles using the 
appropriate Lagrangian Effective Theory (LEFT) [10] and compared these results against numerical 
simulations. Tests of LPT at the density held level have also been performed in [19, 20] and tests of 
Effective Field Theory in connection to LPT have been performed in [21, 22] . Here we use the density 
held induced by those displacements and compare it with the same suite of simulations. We are thus 
testing the original Eulerian Effective Theory of Large Scale Structure [9[. This paper can also be 
seen as and extension of [23] that we use similar techniques and concepts such as the introduction of 
transfer functions. Here the focus is on an interpretation of the results in the EFT framework. 

As we did in BSZ for the displacements, the goal of this study is to test the range of validity 
of perturbative approaches at the held level rather than the level of n-point statistics. We do so 
by comparing iV-body results to the perturbative calculation for the same initial conditions, thus 
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avoiding errors arising from sample variance, i.e. the fact that with the A'^-body code we have only 
simulated one possible realization of the stochastic initial conditions. This is important because some 
of the effects we are after are small and are thus difficult to isolate from n-point statistics in the 
presence of sample variance for reasonable size simulation volumes and numbers of realizations. A 
similar approach was used by [24] to confront SPT with simulations. 

This paper is structured as follows: we will first review the Lagrangian dynamics and the mapping 
to densities in Sec. 2. We will then compute the density power spectra and fields arising from LPT 
and compare them to simulations in Sec. 3. We then extend this framework to estimate its ultimate 
reach by allowing for free transfer functions in Sec. 4. We will uncover an irreducible error that we 
associate with the stochastic term and compare this term to the Lagrangian stochastic term in Sec. 5. 
We will conclude and summarize our findings in Sec. 6. 


2 Dynamics in the EFT framework 


The Eulerian position of particles is given by the sum of their initial Lagrangian position q and the 
subsequent displacement s{q) 

x = q + s{q). (2.1) 

The displacement in LEFT is governed by the equation of motion [10, 18, 21, 22[: 

s + Hs =-V 4>{x) + act + a^toch, ( 2 . 2 ) 

which is solved perturbatively in powers of the density (or its power spectrum) [25[. The gravitational 
potential is related to the density via the Poisson equation Acj) = 3/2f2m'H^(5 and the density is related 
to the displacement via the mass conservation equation [1 + i5(a:)] d^cc = d^q. 

The perturbative expansion is not guaranteed to converge to the correct answer on small scales 
even if all terms in the series are included (see [21[ for an illustration of this in 1 dimension where 
the perturbative theory can be summed). Furthermore, the mistake on small scales affects the large 
scales when statistics at orders higher than tree-level are computed because the momenta in the loops 
can become large. To fix this problem one needs to modify the dynamics. This is accomplished by 
adding the additional acceleration terms and Ostoch on the right hand side of equation (2.2). The 
term is the part of the additional acceleration which can be explicitly computed in terms of the 
perturbative solution of the equations. On the other hand only the statistical properties of Ostoch can 
be computed, and they are not correlated with the deterministic part. 

Symmetries dictate the structure of and agtoch- For example in cosmologies similar to our 
own to compute the one loop power spectrum, the leading correction comes from and has the 
form: act = Z^(t)V(V • s) where I has units of length [10[. It is important to stress that this EFT is 
only meant to be a good description of the dynamics on large scales. On small scales, smaller than 
the non-linear scale, the displacements computed using LEFT are not a good approximation to the 
actual ones. 

In the absence of Oct and Ogtochi in the standard Lagrangian perturbation theory (LPT), the 
scalar part of the displacement field can be solved for recursively in terms of the underlying linear 
density field (i5o): 

^("■^(fc) = do{p^)^6^^\pi+ ...+Pn)ln{Pl,---,Pn) , (2.3) 


where k = pi + ... + Pn and In are kernels than can be found in the literature, e.g. [25-28] . 
The density field arising from the displaced particles reads as 

6{k) = J d^q exp [ifc • (q -|- s(q))] — (27r)^(5^°i (fe) , 

and the density power spectrum is given by 

Ps{k) = 


d^r exp [ifc • r] (exp [ifc • As]) — 1 


(2.4) 
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where As = ^( 52 ) — s(qi) and r — q 2 — qi. The expectation value of the exponential can be evaluated 
using the cumulant expansion theorem 


(exp [—ik ■ As]) = exp — 

+ 


^kikj (AsiAsj) - ^ikikjki (AsiAsjAsi)^ 

^hkjkikm {AsiAsjAsiAsm)c + ■ ■ ■ > 


( 2 . 6 ) 


where the subscript c denotes the connected part of the correlator. 

If one where to expand the exponential in equation (2.6) and keep only terms up to a given 
order in the power spectrum, the resulting series is identical to that of Standard Perturbation Theory 
(SPT)^. In almost all respects these are the only pieces of (2.6) that are trustworthy. The one exception 
are terms associated with large scale motions which are fixed by the equivalence principle and are 
resummed correctly to all orders by (2.6). These terms are of interest in practice as they are the ones 
responsible for the smoothing of the BAO peak. The other terms kept by the exponentiation are a 
superset of the terms one should keep to be consistent at a given order. As long as these terms are 
small, keeping them does no harm and we will mostly do so for computational convenience. If we 
were only to leave exponentiated the bulk motion pieces, one would end up with formulas analog to 
those in the so-called IR-resummation procedure [14, 29]. We will see that in practice keeping terms 
in the exponential is not particularly harmful. 

Even though the trustworthy part of equation (2.6) coincided with what one can compute using 
SPT with the IR-resummation, we found it useful to be able to include the effects of the bulk motions at 
higher orders automatically and at the field level. This is so because in this paper we will determine 
the coefficient in the EFT counter terms by correlating the A-body results with the perturbative 
calculations. If one does this directly in Eulerian space, one is effectively comparing the final density 
with the initial one, measuring so-called propagators [30, 31[, which are significantly affected by the 
bulk flows. The EFT terms are a small correction in these statistics. Performing the comparison 
starting from the Lagrangian displacements we will not have to deal with this complication. 

Equation (2.6) illustrates an important difficulty in going from Lagrangian to Eulerian space. 
Even if one is using the effective theory approach to compute the displacement, the relation between 
displacement and density involves expectation values of quantities evaluated at the same point in 
configuration space. Arbitrarily high momenta contribute to these expectation values and thus the 
Lagrangian to Eulerian transformation introduces new mistakes that need to be fixed by additional 
counter terms. As an illustration let us consider the quadratic term in equation (2.6): 


— -kikj (AsiAsj) 


Ik' ^ - cos(fe' • r)] P{k) 


(2.7) 


Here P{k) is the spectrum of the displacement divergence Si^i. In the Zel’dovich approximation 
the power spectrum coincides with the linear matter power spectrum and in one loop LPT P{k) = 
Plin + 2Pi3^l + T22 ,L) where 2 P 13 l and P 22 ,l are the constituents of the one loop displacement 
divergence power spectra and the subscript “L” stands for Lagrangian. In the Lagrangian EFT this 
power spectrum is corrected by a counterterm ak^Pun due to the fact that the displacement sourced 
by the EFT counterterm in Eq. (2.2) correlates with the linear displacement. This term serves to 
regularize the UV-sensitivity of Pi 3 ,l- 

Equation (2.7) has a zero lag contributions proportional to: 

j(K) 

{si{q)sj{q)) = (s,(0)sj(0)) = ^ . (2.8) 

These are sensitive to the UV part of the spectrum which the EFT does not model correctly. This 
is illustrated in Figure 1. The figure shows that even though on large scales the LPT displacement 
agrees with that computed in simulations it ceases to do so on small scales. It also shows that the 
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^For a review of SPT please see [8]. 








Figure 1. Power spectrum of the displacement divergence —ik s (left panel) and ratio of this power spectrum 
to linear theory (right panel). We show the non-linear displacement divergence from the simulations as well 
as the LPT and EFT terms. The 3LPT displacement overpredicts the true non-linear displacement. This 
agreement can be fixed by adding the leading order counterterm (green dotted) until k — 0.2 /iMpc“^, but 
for larger wavenumbers the counterterm clearly makes the agreement between simulations and theory worse. 


range of agreement can be extended by including the EFT counter term but these new terms only 
improve the agreement of the displacement with simulations on large scales. On small scales they 
could even make the agreement of perturbative and non-linear displacements worse, which will in turn 
degrade the large-scale density through the zero lag contributions. 

Figure 2 illustrates the UV sensitivity in the standard LPT calculation if one where to use 
equation (2.5) as is. For example if s = -|- and the loop integrals are computed using 

a cut-off fcjiiax, the result is very sensitive to that cutoff with the power spectrum differing more and 
more from the simulation results as fcmax is increased. 

The UV sensitivity generated by the mapping from Lagrangian to Eulerian arises from all terms 
in equation (2.6) not just the quadratic piece we just discussed. To understand better all the different 
contributions it is easier to start with the expansion of the density field in terms of the displacement. 
We obtain up to third order: 

S(k) « ik,Si(k) - ^kikj[si * Sj](fe) - ^kikjki[si * Sj * si]{k) , (2.9) 

where [a * b] stands for a convolution in Fourier space. Order by order this expression is equivalent to 
SPT, in particular we have at first order ikis\^\k) = 6^^\k), at second order 


jU) = ik,sf\k) - -hkj[sf^ * , 

f \ 1 w . Ik ■pikp2 

F2[P1,P2) = 7^l2[Pl,P2) + X-2-2“ > 

Z Z P2 

where F 2 is the second order Eulerian coupling kernel [8]. At third order we have 


( 2 . 10 ) 

( 2 . 11 ) 


<5(3) = ikis[^\k) - kikj[s[^^ * s^-^^](fc) - ^kikjki[s[^^ * * s|^^](fc) 


f X 1 , / \ , 1 ; / . fe • Pi + P 2 ) fc • P3 ,1 fe • Pi fc • P 2 fc • P3 

■F3(Pi,P 2,P3 = ^^3 Pl,P2,P3) + ^h[Pl,P2) , , ,, -^ ^- 2 - 

3! 2 Pi+P2r Pi 3! pf 


9 9 

pi pi 


( 2 . 12 ) 

(2.13) 
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However, keeping displacements up to a certain order, we automatically generate higher order terms, 
as we will discuss in more detail below. For the next-to-leading order power spectrum we have^ 

Ps{k) = kikj (si|s^-) +\kikjki ([s^ * Sj]\si) + ^kikjkik^n ([si * Sj]\[si * s^]) - ^kikjkik^ ([s^ * Sj * Si]|sm) 


(2.14) 

Let us consider the terms separately. The power spectrum of the displacement divergence A receives 
contributions from Pn, Pi 3 _l and P 22 ,l- In LEFT, Pi 3 _l comes associated with its counterterm ak^Pu 
[10]. For the bispectrum like term we have^ 


Pi3 = 2ikikjki 




(k) s\^\-k)'^ =4P(/c) J ^l2{k,p)h2{p,k - p)P{p) . (2.16) 


For large internal momentum the above integral scales as 




35 




(2.17) 


which is UV sensitive and thus needs to be regularized by another k'^Pn counter term beyond the 
one introduced in LEFT to regularize Pi 3 ,l- To the extent, that the form of the UV sensitivity is 
the same as the one in Pi 3 ,l, this just changes the coefficient of the k^Pn counterterm. As discussed 
above in Eq. (2.12) [s^^^ * is part of and P 13 will thus be a part of Pi 3 ,e- 
The other contribution arising from the three point correlator yields 


P22 = ikikjki 




= 2 y k - p)h2{p, k - p)P{p)P{\k - p|) . 


For large internal momentum the above integral scales as 


(2.18) 



P\p) 

7^4 


(2.19) 


In the EFT language these UV-sensitivities that are not proportional to the linear power spectrum 
and scale with an external are called stochastic contributions [9, 11]. These terms arise from a 
random reshuffling on small scales that conserves mass and momentum ]12, 32]. 

Yet another stochastic contribution is given by the four point function 


C = kikjkikr, 






) = ^ 2 (p. k - p)P{p)P{\k - p|) . 


For large internal momentum the above integral scales as 


1 


k^ 


pHp) 

ti4 


( 2 . 20 ) 


( 2 . 21 ) 


Finally, let us discuss the cross term between the cubic and linear displacements 


P(p) =-fcV2p(fc) . (2.22) 

^For notational convenience we will often omit the momentum conserving Dirac delta + k') and the normal¬ 

ization (27r)^ in our equations when relating expectation values to the power spectrum. Please keep in mind that we 
always mean (27r)^(5^^^ (fc -h k')P(k) = {S(k)S(k')). 

^ Here we are introducing the convolution vertex 


D = -kikjkikj) 




h2(pi,P2) =1 +Pl ■P2 




(Pl ■ P2)^ 
2 2 
PIP2 


(2.15) 
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This term is proportional to the displacement variance of the field and becomes part of Pi 3 ,e- The 
UV sensitivity requires another counterterm of the form k'^Pn, or equivalently, changes the coefficient 
of the Lagrangian counterterm. 

The Eulerian power spectrum is given by: 

‘2Pi3,e = 2 Pi 3 ,L + Bi3 + D ^ 

P 22 ,E =P 22 ,h + B 22 + C 

in particular, the low-fc limits of the upper line combine to the low-A: limit of T’i 3 ,E) which is usually 
regularized by the Eulerian EFT counterterm c^k'^P [11]. 

As we mentioned before, up to fourth order in the density field the above expressions agree with 
one loop SPT. However, by considering the displacement fields up to a certain order, we automatically 
introduce higher order corrections. One of the most straightforward one of these corrections is that 
in the above terms, all occurrences of cr^ will be replaced by + 2 (t^ 22 (with j^g being 

the displacement dispersion arising from P13.L and likewise for cr^22)- This effectively changes the 
coefficient of the k^P part of the calculation and thus modifies the coefficient of the corresponding 
counterterm. Besides these terms, one generates higher order fields in Eq. (2.9) and also higher order 
contributions to the other correlators discussed above. If these are large, they can not be absorbed 
by a k^P counterterm. 

Finally, we can now return to our earlier comment on the effects of the bulk flows in the Eulerian 
vs. Lagrangian statistics. As an example we can consider trying to determine the EFT parameter 
for the one loop calculation. Our strategy in [18] for the Lagrangian displacement was to examine 
Pi 3 ,h- The bulk flows do not contribute to Pi 3 ,l but they do contribute to Pi 3 ,e as can be seen using 
Eq. (2.22). Thus if we were to apply that same strategy for the density as we did for the displacements, 
we would have to account for the effect of the bulk flows. This seems a bit unnecessary given that those 
bulk flows actually cancel. In an Eulerian calculation they cancel between Pi 3 ,e and P 22 ,e- By working 
with the Lagrangian displacements and keeping them exponentiated we avoid this issue altogether. 
Of course there might be other approaches that one could use, but we found this strategy easy to 
implement especially given that we are comparing the density computed perturbatively with that 
computed using an iV-body simulation, so it is very natural to work with particles and displacements, 
which are in close correspondence to the particles in the A^-body simulation. 

Throughout this paper we are using a suite of A-body simulations discussed in BSZ. There are 
two simulation sizes: the M simulations with a box length of 500 /i“^Mpc and the L simulations with 
a box size of 1500 /i“^Mpc. Besides the possibility to check simulation convergence, the different 
size of the simulations also imposes different Nyquist frequencies and thus different cutoffs in the 
perturbative displacement calculations on the simulation grid. We can use this difference to test the 
convergence of our results. Both simulations calculate the gravitational evolution of 1024^ particles 
from the initial redshift z = 99 to present time. We presented numerous numerical tests of our runs 
in BSZ so we will not repeat them here. 

3 One loop EFT calculation 

The mapping from Lagrangian to Eulerian space is a non-linear transformation that is not additive, 
i.e., the effect of the n-th order displacements can not be calculated in isolation, but comes always 
associated with all lower order contributions. Throughout this paper we will thus consider density 
fields obtained by using the displacements computed up to a certain order of LPT and denote the 
highest occurring order by an overbar. Hence, is the density field calculated using all displace¬ 
ments up to n-th order. At first order there is no ambiguity and we will thus often omit the overbar 
for the fields arising from the Zel’dovich displacements. A one loop calculation requires including the 
terms up to third order, so we start by considering^ 

5 = 5 ( 3 )-H (3,1) 

■^Note that the EFT parameter a employed here is should not be confused with the Lagrangian EFT parameter 
employed in BSZ. 
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Figure 2. Performance of LPT density power spectra for ILPT/Zel’dovich and 3LPT evaluated for the same 
initial conditions as the simulation for redshifts z = 0 (left panel) and 2=1 (right panel). We clearly see 
the missing power in both the ILPT and 3LPT power spectra. This failure can be reduced by setting the 
maximum wavenumber or LPT cutoff to fcmax = 0.4 hMpc“^. The same effect can be achieved by adding a 
k^Pii counter term as predicted in the EFT framework. The jaggedness of the ILPT lines compared to the 
higher orders of LPT arises is due to the missing higher order terms that also cancel cosmic variance. 


This displacement contains terms that go beyond one loop Eulerian EFT. Because the displacement 
has been kept in the exponent, the expression contains high powers of and that are of 

the same size as terms being dropped, for example from When computing power spectra using 
<5^^) we also include the correlator that is formally a two loop term. The fact that we are 

keeping higher order terms is no particular problem as we are not expecting Eq. (3.1) to be correct 
at higher order, and will only trust it in the regime where these extra terms are negligible. 

To obtain equation Eq. (3.1) we have also kept only terms linear in the counter terms. This 
approximation should be correct in the same sense, in any case we cannot trust higher order terms 
arising from the exponentiation of the counter terms as higher order counter terms contribute with 
similar magnitude. 

One potential source of worry is that the counter terms are needed to fix mistakes in the pertur¬ 
bative calculation which are kept in the exponential, while the counter terms themselves have been 
brought down. This would be a problem if there was a very large cancellation between the pertur¬ 
bative terms and the counter terms that will no longer happen after the approximation. Fortunately 
this is not the case at one loop, the counter terms only make a small difference in that case and the 
inaccuracy from the expansion is not larger than the two loop terms we are neglecting at this order. 

Using Eq. (3.1) we can calculate the auto power spectrum of the model and the cross power 
spectrum between the model and the linear (Zel’dovich) density field, 

P^i i = Pi 3 -I- ak^Pii and Pni = ^33 + 2ak^Pi^ + a^k^Pn . (3.2) 

To the extend that the a^k^Pn term is a higher order correction, we will drop it in the following. 
The cross correlation has the advantage, that we can explicitly probe the phases of the counterterm. 
This might help to avoid overfitting that could happen when using only the final power spectrum to 
measure the coefficients, where they could be degenerate with higher order corrections. 

For reference, we will also consider the Zel’dovich density and associated power spectrum: 

P„i j = {l + ak^)Pii , and P„i = {1 + ak^fPn . (3.3) 

Note that we have included a counter term for Zel’dovich. As we discussed before, this is needed 
because even though Zel’dovich is a linear calculation, the Lagrangian to Eulerian mapping introduces 
zero-lag terms (for instance through Eq. 2.22) that need correction. 

We show the performance of the Lagrangian EFT models in Fig. 2. We clearly see that while 
the bare LPT power spectra reproduce the full non-linear power spectrum on the very larges scales. 
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Figure 3. Scale dependence of the coefficient a of the leading order counterterm in the L simulation (left 
panel) and the M simulation (right panel). We consider the one loop LPT model (green, lower) and the 
Zel’dovich model (red, upper). The solid lines show constraints from the cross spectrum with the Zel’dovich 
field, whereas dashed lines show constraints from the auto spectrum. The horizontal gray lines indicate the 
value of the EFT coefficient employed in Fig. 2. The difference between the green and red lines on large 
scales arises from the low-fc limit of Pi 3 ,l (dS/BSirJfc^P) as well as and 0 ^ 22 entering in Eq. (2.22). The 
difference between the green lines in the two panels on large scales is given by the cutoff dependence of crj 13 
and cr|^ 22 - 


they lack power on smaller scales. Despite the lack of power, the ratio of LPT and simulations is 
smooth, which means that IR motions have been appropriately resummed and that the BAO wiggles 
are well reproduced. Note in particular that the plots show the results of a single realization, i.e., by 
performing the LPT calculation for the same seeds we cancelled cosmic variance (see [24] for a similar 
technique employed to test SPT). This cancellation of higher order cosmic variance contributions is 
apparent when going from ILPT to higher order LPT, for which the scatter is considerably reduced. 
This is due in particular to odd correlators (for instance P 12 ) that vanish only when averaging over 
many realizations. 

Fig. 2 also shows that the ratios are highly cutoff dependent, especially for 3LPT with higher 
cutoffs leading to a stronger suppression due to the artificially large small scale displacement power in 
3LPT. In this context, it is illustrative to look at the contributions to the r.m.s. displacements from 
modes of different scales, which are shown in Fig. 4. In turn, by choosing an appropriately low cutoff 
one can bring the 3LPT prediction and the non-linear power into agreement up to /c = 0.1 hMpc”'^. 
Let us now consider the EFT corrections to 3LPT, by choosing a in Eq. (3.2) appropriately we can 
improve the agreement between theory and simulations on large scales and obtain a one percent fit up 
to fc = 0.1 /iMpc”^. For comparison we also consider an EFT counter term on the Zel’dovich power 
spectrum only, which performs almost as well. 

The measurements of a as a function of scale is shown in Fig. 3 for the L and M simulation. Let 
us first consider the counterterm on the Zel’dovich power spectrum shown by the upper pair of red 
lines. We consider the constraints from the cross-power of Zel’dovich and the non-linear density field 
and the constraint from the auto power spectrum of the density 

_ ^nl,l ~ Pii _ Pni — Pii 

across — ’ aauto — 2k'^ 

which are depicted by red dashed and solid lines respectively. At the level of the displacement field 
they would correspond to the error and non-linear estimators employed in BSZ. We see that for both 
the L and M simulations the a parameters from both the auto and the cross estimator quickly decay 
from their initial amplitude of a « 6.5 h“^Mpc^. This scale dependence or “running" of the EFT 
parameter is a first indication of the breakdown of the ansatz and the presence of higher order terms. 
A further indication of higher order terms playing a role is the fact that the auto and cross estimators 
deviate. Finally, let us comment on the size of a we are finding and contrast it with the typical 
















Figure 4. Cumulative contributions to the displacement dispersion up to a certain maximum wavenumber 
femax from the various contributions to one loop LPT. We show both the bare contributions of LPT (solid) 
as well as the ones from the displacements regularized by the transfer functions (dashed). The non-linear 
displacement dispersion coincides with the linear one before a transfer function is employed. The transfer 
function has a mild influence on the 22 '’’d.ii part, but significantly reduces the contribution from 5 ^ 3 . 


value Of® « 1 usually found when comparing the EFT prediction with the power spectrum 

of the density computed directly in SPT. In Zel’dovich the linear power spectrum is damped with a 
factor k^a^/2 = 18k^ /i“^Mpc^, whereas in SPT one gets 61/210 /c^ct^ « 10.45fc^ ft,“^Mpc^. The true 
damping is roughly (61/210tT^-|-a®)fc^ « 11.45fc^ /i“^Mpc^ and thus one expects a = 6.5 /i“^Mpc^ for 
the Zel’dovich transfer function, which is what we see in Figure 3. This result is fairly independent 
of the cutoff scale since the Zel’dovich displacement dispersion converges quickly (see for instance 
Fig. 4). 

Let us now discuss the 3LPT case, for which the estimators read 


Q^cross 


^nl,l ~ -P13 

k^P-ri 


<^auto 


^nl ~ P33 
2BPll 


(3.5) 


They are shown as the green lower set of lines in Fig. 3, and show a flattening on large scales,® which 
we would indeed expect over the range of validity of the ansatz in Eq. (3.2). For both simulation sizes 
and cutoffs the constraints start to deviate from a constant and the cross and auto constraints from 
each other &t k k. 0.08 /iMpc”^. This scale actually coincides with the scale where we expect two 
loop corrections in Eulerian Perturbation Theory to matter [33]. The difference between the low-fc 
limits of the coefficients in the L and M simulation can be explained by the strong cutoff dependence 
of cr^ 22 and ^ 3 . 

In the right panel of Fig. 2 we show the performance of the LPT and EFT terms for redshift 
z = 1 , for which our fiducial cutoff of /cmax = 0.6 /iMpc”^ leads to almost optimal results, providing 
a one percent accurate fit up to A: = 0.2 /iMpc”^ without a counterterm. As we saw from the z = 0 
example, a similar result would have been obtained with a different explicit cutoff and an appropriate 
counterterm. Thus there is nothing special about this particular cutoff. All we are seeing is the 
familiar phenomenon that the size of the counterterms runs with the cut-off and redshift. In fact the 
counterterms are what allows the theory to make predictions that are cut-off independent. 


4 Going beyond one loop 

In BSZ we compared the displacements calculated with LPT to those produced by an A^-body code. 
Although on large scales there was excellent correlation, they differed significantly on small scales. In 

®The precise extraction of the EFT counterterm coefficients on large scales is complicated. While being the cleanest 
place to extract a low energy constant, the counterterm is only a tiny correction 0(10“^) a,t k = 0.03 /iMpc“^, and 
simulation codes are struggling to reproduce linear growth at the 0(10“^) level [6]. 
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a similar fashion as [23] we defined transfer functions ai{k) such that 

spT(fc) = ai(fc)s^^^(fc) + a2{k)s^‘^\k) + a^(k)s^^\k) H- (4.1) 

effectively minimizing the difference between the fV-body answer and spx- These transfer functions 
describe the motions induced by a^t in LEFT as well as higher order contributions. Their scale 
dependence is shown in Fig. 18 of BSZ. The full 7V-body displacements are given by 

body ~ ^PT T ^stoch; (4-^) 

where Sgtoch results from Ogtoch in LEFT. 

The transfer functions on LPT displacements beyond the leading order go to zero on small scales, 
indicating that the LPT displacements are wrong on these scales. Without these transfer functions 
one-loop LPT (in particular its Pi 3 ,l part) significantly overpredicts the r.m.s. displacements. But 
also the Pn.p r.m.s. displacement is overpredicted by 1.7 /i“^Mpc^, and it is a mere coincidence, that 
the bare value (without the appropriate displacement transfer function) agrees so well with the true 
non-linear r.m.s displacement. Fig. 4 shows the different contributions to the r.m.s. displacement. In 
BSZ we were also able to show that with the transfer functions we include LPT contributions of orders 
exceeding the one explicitly computed. For example both and contributions are captured 

by including the Oi transfer function. We showed that 3LPT with transfer functions (3tLPT) contains 
all the terms included in the two loop LPT calculation. Furthermore, we saw that fourth order LPT 
didn’t improve the agreement with simulations because at that point the main source of error came 
from Sstoch- With the transfer functions we consider in this paper we are capturing all the terms in 
the two-loop LPT calculation. 

In this section we will consider the density fields generated from displacement fields with transfer 
functions. As we have argued above, the presence of additional zero lag terms in the mapping between 
Lagrangian and Eulerian space introduces additional mistakes. The symmetry of the EFT constrains 
the structure of the additional terms needed to fix those mistakes. In this section we will include 
these terms by introducing an additional transfer functions to the density field. The fact that the 
density field computed from the best possible displacements differs from the full results in ways that 
can be absorbed by a simple transfer function was already noted in [23]. Here we extend the LPT 
calculation to higher orders and also make the connection with the Eulerian EFT calculations. 

As in Lagrangian space, as a metric of the performance of a model for the density field (rather 
than the power spectrum), we consider the error power spectrum defined as 

Terror ~ (l^nl ^model I ) ~ (^nl|*^nl) “f (^model l^model) 2 (<5^1 [(^model) 7 (4-^) 


where here model refers to the perturbation theory calculation computed up to a given order. On a 
given density model we can define an overall density transfer function by 


T = 


(^nll^model) 
(^model I ‘^model) 


(4.4) 


such that the optimal error power spectrum and the ratio of error and non-linear power spectrum are 
given by 


Terror, TF — (<^nl|^nl) 


(^iill^model) 

( ^model I ^model) 


(4.5) 


Terror,TF .. (*^nl [f^inodel) _ -1 2 f ^ 

Pnl ~ ~ (^„l|<5nl) (5model|<5model) " " ^ ^ 

where Tcc is the cross correlation coefficient. This statistic measures the r.m.s. deviations between 
a certain model and the non-linear data. It is quantifying the performance of the model at the field 
level rather than the level of the power spectrum. 

Fig. 5 shows the error power spectrum for different perturbation theory calculations. As an 
illustration we show the results of LPT without any transfer functions (EFT counterterms). One 
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Figure 5. Error power spectra at redshift z = Q for the L simulation and various orders of LPT and tLPT. 
Left panel: Total power. The gray lines show the expectation for the stochastic term for fcni = 0.3 /iMpc“^ 
and fcni = 0.5 /iMpc“^, respectively. We see a decrease in mode coupling as we go to higher orders and 
implement transfer functions on the density. Right panel: Ratio of the error power spectrum and the full 
non-linear power spectrum. 


striking observation in this case is that including makes things worse. The curve labeled T5{\ + 2) 
is effectively the equivalent of the LEFT one loop calculation. Even if is not included, the only 
part of relevant at this order is the cross term 1 — 3 which is captured by the transfer function. 
We have also computed another version of an effective one loop calculation T(5(l -h 2 -f 3) which gives 
very similar results (although the transfer function is different). The comparison between TS(1 + 2) 
and 5(1 -f 2 -f 3) illustrates the fact that the overall transfer function is indeed crucial, it is better to 
include that than to add an additional order in the displacement field if the goal is to minimize the 
error at the level of the density field. 

The line labeled TS(T1 + T2 + T3) uses transfer functions for all the displacements up to the 
third order and plus an overall transfer function on the density to correct the leading order part of 
the zero-lag terms from the mapping. This model yields the lowest error power spectrum, such that 
we interpret it as the best possible perturbative model. We have also compared this to the case 
where we have several transfer functions at the level of the density determined in the same way as the 
displacement transfer functions in BSZ: T35(Tl-|-T2-|-T3)-|-T25(Tl-|-T2)-|-ri5(n) (see App. A). The 
latter contains all the terms in a two-loop EFT calculation and gives very similar results to the cases 
considered here.® The error we see in T5(T1 -|- T2 -|- T3) does not decrease appreciably if we include 
higher order displacements either. We interpret this error as arising primarily from the stochastic 
displacements Sstoch- 

We had uncovered a stochastic contribution to the displacement already in BSZ. The stochastic 
term in the EFT arises from a mass and momentum conserving shuffling of mass on small scales. As 
a result, its power spectrum has to scale as for low wavenumbers and we saw this behavior clearly 
in BSZ for the divergence of the displacement. Fig. 5 however indicates that this behavior seems to 
be violated by the error power spectrum in our simulations even on rather large scales, which we are 
identifying with the stochastic term. We have performed a number of numerical tests to check the 
stability of this stochastic contribution and found it to be stable. We show some of the checks in the 
Appendix, for example Fig. 10 shows the comparison of the error power spectra found when analyzing 
the L and M simulations. 

We will discuss the shape of the error power spectrum in more detail in the next section. What 
we will see is that on very large scales, this stochastic term of the density indeed agrees with the power 

®At the two loop level, the linear counter term enters correlated with itself as Pjj, but the second and third order 
counterterms enter only correlated with perturbative terms as ^22 and To this extent, the counter term if present 
in the data, should correlate with the perturbative second and third order basis vectors. Thus our T3(5(T1 + T2 + T3) + 
T2(5(T1 + T2) + TiS(Tl) captures all relevant terms for two loops. In fact one expects this example to be if anything 
better than one might do in a first principle Eulerian EFT calculation as we are allowing the transfer functions to have 
an arbitrary shape. We will discuss this in more detail in App. B. 
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Figure 6. Left panel: Ratio of the best possible EFT power spectrum to the non-linear power spectrum as a 
function of redshift. We indicate the 1% and 10% accuracy lines and mark the crossing of the 1%-threshold 
by vertical lines, whose wavenumbers are given in Tab. 1. Right panel: Ratio of the perturbative model with 
and without transfer functions and the non-linear power spectrum at 2 = 0. 


spectrum of the stochastic displacement divergence, as it should. But we will also identify corrections 
that arise from the mapping from Lagrangian to Eulerian space leading to deviations from the 
scaling on surprisingly large scales {k « 0.03 /iMpc”^ at z = 0). 

In Fig. 6 we show the ratio of error and non-linear power spectrum Pen/^NL for three redshifts 
z = 0,1, 2 for the TS{T1 + T2 + TS) example to quantify up to which wavenumber the perturbative 
calculation can be expected to agree with the A^-body result. We quote the wavenumbers at which 
the stochastic power crosses the 1% and 10% level in Tab. 1. While one should not focus too much 
on the specific values, one should definitely note the steepness of the curves in the left panel of Fig. 6. 
This means that at a fixed k away from the non-linear scale, the size of the error changes dramatically 
as one goes to higher wavenumbers. This is important, since for data analysis applications, such as 
trying to see the small effects of primordial non-Gaussianity in the two- and three point functions [2] , 
precision will probably be more important than reach. 

The right panel of Fig. 6 shows Pmodei/PNL — 1 and illustrates again that the biggest improvement 
in reach comes from the inclusion of the final transfer function, fixing the problems caused by the 
mapping. The comparison between T5{\-\-2) and r(5(l-|-2-|-3), which are both effectively equivalent to 
one-loop EFT calculations (with higher derivative counterterms) shows the difference that the higher 
order terms that are only partially included can make. It is amusing to note that T(5(l -I- 2 -|- 3) is 
actually slightly worse, so the additional work to include did not result in an improvement here. 
This is perhaps not surprising given how bad is on small scales and the fact that the Lagrangian to 
Eulerian mapping makes the large scale density depend on these mistakes. Of course with additional 
freedom from more counter terms one should be able to absorb these differences. Finally one may 
notice that in terms of reach, T5{T1 -\- T2 + T3) does not even improve by a factor of two. But reach 
is perhaps the wrong metric as the error curves are very steep. Fig. 5 shows that away from the 
non-linear scale, the error in T6{T1 -|- T2 + T3) is smaller than the one in T5{1 + 2) by about one 
order of magnitude. 

A map of the various density fields discussed in this section is shown in Fig. 7. It clearly shows how 


z 

h% 

^10% 

0 

0.25 /iMpc"^ 

0.46 /iMpc”^ 

1 

0.48 /iMpc"^ 

0.98 /iMpc”^ 

2 

0.85 /iMpc”^ 

1.72 /iMpc"^ 


Table 1. Wavenumbers, where the stochastic term amounts to a 1% or 10% correction to the non-linear 
matter power spectrum. 
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Figure 7. Non linear transformation of the density field in a patch of 300 h ^Mpc length and 15 h ^Mpc 
depth. 


well correlated the structure in a Zel’dovich realization is with the non linear structure. Overdensities 
are washed out and voids are clearly underdense. Adding higher order displacement fields and transfer 
functions on the displacement fields has no strong imprint in this picture beyond a slight sharpening 
of the overdensities and filaments. The final density transfer function shown in the last panel clearly 
has the strongest effect, most remarkably a sharpening of structures in all environments. But even at 
this level there are still obvious differences between the best perturbative approach and the non-linear 
field. 


5 Stochastic Term 


In BSZ we identified an irreducible error at the field level that we associated with the stochastic 
term of the EFT. In this Section, we are relating the Lagrangian stochastic term to the Eulerian one. 
From now on we will consider displacement fields up to a certain order including all transfer functions 
and denote them Spx, in particular we will be mostly concerned with the displacement fields up to 
third order, he., Spx = -I- The total displacement field is then the sum of the 

perturbative and the stochastic part s = spx + Sstoch- 


Toy model: origin of the transfer function for the density 

Let us consider the case where we expand only the stochastic displacement in Eq. (2.4) 


{2TT)^S^^\k) + 6{k) « y d^gexp [ik ■ {q + spx)] + ifciSstoch.i - ^hkjSstoch,^Sstoch,: 


(5.1) 


21 ^stoch,2'^stoch,j ^stoch,/ 
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The stochastic term is likely related to small scale non-linear phenomena and is orthogonal to the 
perturbative terms by construction. Averaging the above equation over the short modes and assuming 
that there are no correlations between the statistical properties of the stochastic term and the long 
modes we get: 

(27r)^(5(°)(fc) -k « J d^gexp [ifc • (q -k spT(g))] (^1 “ ^ ((^ ' (5'2) 

This result could have been obtained at the level of the exponential using the cummulant expansion 
theorem and neglecting coupling between short and long modes 


(27r)35(D)(fc) + (5(fc)),,„^, 


d^gexp [ifc • (q + Spx)] exp 


1 

2 


k'^a 


2 

djStoch 


(5.3) 


and thus 


6{k) = 6pT{k) exp 


1 

2 


k^aj 


,stoch 


(5.4) 


This is a Gaussian smoothing of the PT results, it has the form of a quadratic counterterm and its 
amplitude is given by the stochastic displacement cr'^stoch ~ /i“^Mpc. It also indicates, that for 
the best density field one should consider the best possible result for the displacement field and then 
multiply with an overall transfer function that captures the r.m.s. displacement of the stochastic term. 
Of course there are additional contributions when one takes into account that the statistical properties 
of the stochastic term depend on the long modes. At the lowest order this was discussed already in 
[10] and found empirically in [23], here our goal was to illustrate this point again for convenience of 
the reader. 


Mapping of the stochastic terms 

Let us reconsider the mapping between Lagrangian displacements and Eulerian density Eq. (2.4), 
distinguishing between PT and stochastic terms 


{2T:f S^^^k) + 5{k) K J d^gexp[ifc-q] ( 1 + ifci(spT + Sstoch)* - ^fc*fci(spT + Sstoch)? 


1 


ikikjki{spT 


'5stoch)ij/J ■ 


(5.5) 


The relevant contributions correcting the Lagrangian stochastic term are then given by 

Pstoch,E(^) =kikj (Sstoch.il^stochj) 2\kikjki ([Sstoch.i * l-^stoch,/) ikikjkl ([5stoch,i * -^stoGhj] l-^stoch,/) 

+kikjkikm ([Sstoch.i * SPTj]|[Sstoch.i * SPT.m]) + hkjklkm ([spT.i * SPTj * Sgtoch,/] |Sstoch.m) 

(5.6) 


Using the fact that the mode coupling is orthogonal to perturbation theory, we can readily calculate 
the contributions from the even correlators in the above calculation. In particular, we have for the 
first, fourth and fifth term 


Pstoch,E(^) UPstoch,L(^) “t” Pstoch,L(^) 




k ■ p k ■ {k — p) 


[k — pY 

=Pstoch,L(fc)(l - + / ^2(P, fc-p)Pstoch,L(p)Pll(fc-p), 

J -P 


Pstoch,L{p)Pll{k - p) 


(5.7) 


where we assumed Gaussianity for the stochastic term. Unfortunatly, this calculable correction is 
too small to explain the differences between the two stochastic terms identified in the simulations in 
Lagrangian and Eulerian space. 
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Figure 8. Contributions to the mapping from Lagrangian to Eulerian space. We start from the Lagrangian 
mode coupling (green dash-dotted), to which we add the next to leading corrections (Satochl[Sstoch * Spt]) and 
(Sstochl[Sstoch * Sstoch])- There is a significant cancellation between these two terms themselves and with the 
original mode coupling term. The sum of these three terms (black) traces the Eulerian mode coupling defined 
in the previous section (red dot-dashed) very well on large scales. 


In BSZ we were only concerned with the two point functions involving the fields and not the three 
point functions. Thus, we only established the orthogonality of stochastic term and perturbative terms 
at the level of the two point function and did not consider the three point function of the stochastic 
field itself or possible three point functions between the stochastic and the perturbative fields. These 
could arise if for example, the stochastic piece is modulated by the perturbative modes, leading to a 
coupling ' 

= / K^ip, k - p)4och(p)<^PT(fe - p) ■ (5.8) 

Jp 

with an unknown coupling kernel K 2 . This term does not correlate with pure perturbation theory 
terms and does thus not violate our requirement that (^pxl^stoch) = 0. Alternatively there can be 
a correlation of the stochastic term with the square of the stochastic term, that can be encoded by 
another kernel and the replacement 5px —>■ ^stoch in th® above equation. This is indeed quite natural 
for a non-Gaussian field generated by non-linearities. With the above couplings we would have for 
the three point functions 

(<5pT<5stoch|4toch) = / K 2 {p,k - p)Pstoch{p)Pii{k - p) + P^tochik) / K 2 {k,p)Pii{p) (5.9) 
J p J p 

We have no a priory knowledge about the coupling kernels or the three point functions of the stochastic 
term and thus have to extract the three point function from simulations. 

We measured the correlators ((5pTi5stoch|<5stoch) and ((5stock<5stoch|<^stoch) for the stochastic term 
identified in our simulations and show the results in Fig. 8. We clearly see that the Lagrangian stochas¬ 
tic term for the displacement dispersion differs from the Eulerian one on all but the largest scales. 
In particular, the Eulerian stochastic term is smaller than the Lagrangian one for k < 0.25 ft,Mpc~^ 
and exceeds it for higher wavenumbers. This behaviour is qualitatively expected from the mapping 
between Lagrangian and Eulerian coordinates and the collapse of structure. More quantitatively, we 
consider the corrections from the correlator of the stochastic field with the product of perturbative 
and stochastic fields at a different location and the correlator of the stochastic field and the square 
of the stochastic field at another position. The first of these two terms, leads to order unity negative 
corrections to the original stochastic power around k « 0.1 /iMpc”^, such that only the bispectrum 
of the stochastic term survives. The sum of the original term and the two leading corrections yields 
a result that agrees with the stochastic term identified in the previous section, by projecting out the 
perturbative part of the density field. This is on the one hand very reassuring, since we have recovered 

^For illustrative purposes we will use the shorthand notation <5 = —ifc ■ s. 
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the stochastic part of the density field in two independent ways. On the other hand, the importance 
of three point functions in the mapping from Eulerian and Lagrangian makes it very hard to assess 
the size of the Eulerian stochastic term based on the Lagrangian two point function only. 

The most important observation from the EFT point of view is that the stochastic term asymp¬ 
totes to its expected leading scaling only on very large scales k « 0.01 hMpc“^ and has a shallower 
slope where it matters at the percent level k « 0.2 /iMpc”^. This complicates its modelling in the 
EFT framework, since higher orders in derivatives or additional scales besides the non-linear scale 
need to be considered. 


The stochastic term in the halo model perspective 

The halo model (for a review see [34]) splits the mass distribution in the Universe into distinct 
subsets or regions, each with its center of mass and a density profile. The correlations are then 
disentangled into correlations between these objects and correlations within the objects. Since the 
objects are often taken to be dark matter haloes, these terms are known as one- and two-halo terms, 
respectively. Clearly, the correlations within haloes are non-perturbative and the corresponding halo 
profiles are thus usually extracted from numerical simulations or fitting functions to the latter. The 
only perturbative part is the correlation between the centers of the overdense regions. In perturbation 
theory, there is no such distinction. On the contrary, whenever one employs models for the clustering 
of the dark matter field, one predicts the field at all positions, even within highly overdense regions, 
such as dark matter haloes. Thus the perturbation theory result contains a one halo term, that is 
clearly wrong, since there is no hope that a perturbative approach will recover the virialized structure. 
This was previously discussed in detail in [35] . Here we would like to discuss the sizes of the terms in 
the context of our measurements. 

The power spectrum of terms within a single halo is given by 

r A/r2 

-PiH(fc) = j dMn{M)—u^{k) , (5.10) 


where n{M) is the halo mass function and u is the normalized profile satisfying u{k) -1. The 

amplitude of the one halo term is thus given by the mean squared density fluctuations 


Hih = J dMn{m)^^ « 330/i ^Mpc^ . 


(5.11) 


The above integral is dominated by massive dark matter haloes M « 10^"^ h~^MQ. 

Assuming that perturbation theory correctly predicts the center of mass of haloes, the remaining 
profile error between the perturbative and non-linear profiles upx and unl will contribute the following 
mistake 


f r 12 

APm = / dM n{M)^- wnl]^) - upT(fc) 


(5.12) 


Imposing mass and momentum conservation on unl ~ wpT, the kP and k^ components cancel from 
this term and it thus starts as Thus we have finally that APih oc k'^ as one should expect for a 
mass and momentum conserving local process. This does not mean that the one halo term by itself 
is compensated, just the difference between the PT and non-linear profiles is compensated. 

We identified haloes in our L simulation using a Friends-of-Friends halo finder with linking 
length 0.2. We then extracted the non-linear halo profiles and the profiles of the particles displaced 
by Zel’dovich rather than non-linear dynamics. Phenomenologically, the halo constituent particles are 
roughly at the right position but the PT halo is more dispersed. Fig. 9 shows the profile differences 
between the true non-linear profiles and the Zel’dovich profiles as well as the non-linear profile itself. 
Above fc « 1 /iMpc”^ the two agree, but on larger scales the compensation kicks in and the the profile 
difference asymptotes to the expected k'^ behaviour. If modeled according to the self similar scaling 
Astoch = {k/kn\f [13], this term corresponds to fcni = 0.5 hMpc”^. This profile difference is clearly 
a lower limit on the size of the stochastic density power spectrum since other non-linear structures 
(walls, filaments and voids) and the stochastic error on the center of mass will contribute as well. 
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Figure 9. Power spectrum of the non-linear one halo term Eq. (5.10) extracted from the simulations and 
the one halo term of the profile difference between the simulation and ILPT Eq. (5.12). The gray line shows 
1% of the linear power spectrum. We can clearly see that the profile differences amount to percent level 
corrections to the linear power spectrum at fc « 0.3 /iMpc“'^. We also show the stochastic term in Lagrangian 
space and Eulerian space. The Eulerian stochastic term is larger than the profile difference on large scales 
but approaches it at fc « 1 hMpc“^. 


Given the large scales/small wavenumbers where the stochastic term asymptotes to it is clearly 
not dominated by dark matter halo profiles at these scales. Yet, the scale where the halo profile 
difference amounts to one percent of the linear power spectrum is k = 0.28 hMpc”^, which is close 
to the scale where the stochastic term crosses the 1% threshold. Recently, [36] found a compensated 
one halo term that is orders of magnitude larger than the term discussed here. This difference arises 
from the fact that we are taking the difference between ILPT and the simulations at the field level, 
whereas it was taken at the level of power spectra in [36]. 

6 Conclusions 

Using the IR resummed basis provided by Lagrangian Perturbation Theory, we compared the pertur¬ 
bative density fields with the results of A^-body simulations sharing the same initial conditions. 

Using LPT generated displacement fields with a transfer function on the linear field, as 
suggested by the EFT at leading order, we manage to model the non-linear power spectrum to 1% 
precision up to fc « 0.1 /iMpc”^ at z = 0. Precise measurements of the EFT coefficient show a scale 
dependence of the coefficient extracted at 3LPT level for k > 0.07 /iMpc”^ that is likely due to the 
presence of missing two loop corrections. We also showed that the leading EFT counterterm is able 
to capture the cutoff dependence of the LPT density field for the two cutoffs fcmax = 0.6 ft,Mpc~^ 
and fcmax = 2.4 /iMpc”^ considered here. To assess the maximum range of validity of perturbative 
approaches for the density field and in order to avoid issues with small scale spurious motions we then 
employed the regularized displacement fields with transfer functions defined in [18] . 

We find that in our highest order calculation, at redshift z = 0 the power spectrum of the density 
field is reproduced with an accuracy of 1% (10%) up to A: = 0.25 hMpc”^ {k = 0.46 /iMpc”^). We 
believe that the dominant source of the remaining error is the stochastic contribution, orthogonal to 
the perturbative basis. The stochastic term will likely put an upper limit on the range over which the 
non-linear power spectrum can be modeled by IR-resummed two-loop Eulerian EFT. Perturbation 
theory approaches should generally only aim to model the deterministic (non-stochastic) part of the 
density field. 

The Eulerian stochastic term deviates from the Lagrangian stochastic term of the displacement 
divergence on all but the largest scales. This is expected from the collapse of structure and can 
be explained quantitatively by the three point function of the stochastic term itself and correlations 
between the square of the stochastic term and perturbation theory. The stochastic term only scales 
as on the largest scales, being shallower in the range where it starts to affect the power spectrum 
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at percent level and finally asymptoting to the halo profile. These deviations from the simple k'^ be¬ 
haviour over the relevant scales will greatly complicate the modelling of the stochastic term, requiring 
additional parameters or scales besides the non-linear scale. 
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A Cross checks 

In the left panel of Fig. 10 we compare the error power spectra between the L and M runs. Despite 

the varying simulation resolution and the difference in the cutoffs employed for the perturbative 

calculation the results are very well converged on the scales of interest. In particular, the scale where 

the stochastic term corresponds to one percent of the non-linear power spectrum agrees between the 
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Figure 10. Left panel: Comparison between ratio of error and non-linear power spectra for the M and L 
simulations before and after employing the transfer functions. Despite the very different box size and cutoff 
wavenumbers of the LPT calculation, we find extremely good agreement between the error power spectra. 
Right panel: Comparison of the error power spectra in the M simulation with three transfer functions on the 
displacement fields and an overall transfer function on the density field with the cases where we add the fourth 
order displacement with the corresponding transfer function and the case where the density fields generated 
from the first, second and third order displacements have different density transfer functions. 


two cases. We can thus conclude that we have employed sufficiently general transfer functions to 
capture the cutoff dependence of the theory. 

In the right panel of Fig. 10 we show the error power spectrum for a few cases not explicitly 
considered in the main text in order not to confuse the reader. We start from our fiducial model 
T5{T1 -\-T2 + T3) and generalize it with additional transfer functions. As a first approach we add the 
fourth order displacement field in Lagrangian space multiplied with its transfer function. Considering 
the error power spectrum arising from the corresponding density field with an additional density 
transfer function, we see no considerable improvement over the fiducial case. As the most general 
case we split the basis for the density field, considering the basis vectors <5(T1), S{T1 + T2) and 
5{T1 -b T2 -b T3) separately and allowing three separate transfer functions for these density fields. 
This model has in total six free functions of wavenumber (three in Lagrangian and three in Eulerian 
space). Even in this fairly general approach, the improvement in the A:-reach is negligible. However, 
one very large scales the error at a fixed k decreases by about a factor of two. For the latter term we 
can explicitly show, that it should contain all the terms in two loop SPT and their counterterms. The 
1-5 contribution is accounted for by a transfer function on i5(Tl), the 2-4 contribution is accounted 
for by a transfer function on the 6{T1 + T2) — 6{T1) and the 3-3 contribution is accounted for by 
S{T1 + T2 + T3) — (5(T1 -b T2) for which we actually wouldn’t even need a transfer function to get to 
two loops. Concerning the counterterms, they come in the form of I — 1, 2 — 2, 3—1 and 1 — 1 (where 
the tilde stands for a counterterm involving density fields of this order). All of these terms have been 
included by the transfer functions that we are allowing for. 

B Equivalence with 2-loop SPT 

When expanded SPT and Lagrangian solutions are equivalent. This means that 

J(5) ^ J(l) J(2) ^(3) j(4) ^(5) . (3^^) 

where the dots stand for terms that are higher order than what we need at two loops. We can now 
rewrite: 


j(5) ^ (^(5) _ ^(4)^ ^j(4) _ ^(3)^ (-j(3) _ j(2)^ ^j(2) _ ^(1)^ ^(1)^ 

= -b + (5(3) -b (5(^) -b (5(3) 


(B.2) 
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with the definition 5*-"^ = 5*-"^ — We note that starts at order n, but includes higher orders. 

In fact it correctly includes the effects of large scale bulk motions and this is why it is the convenient 
basis to compare directly with the density of Wbody simulations by means of cross correlations to 
determine transfer functions (or EFT counter terms). 

Because of the equivalence between SPT and the Lagrangian treatment when expanded we also 

have: 

P55 = P2-I00P, E + ■ ■ • (B. 3 ) 

where again the dots stand for higher order terms. Now we can use equation (B.2) to write: 

Plb = Pii + P 22 + ^33 + ‘^Pi2 + 2^13 + 2Pj 4 + Pjg + 2 P 23 + 2^24 + ■ ■ ■ (B.4) 

where we have left out any higher order terms, such as ^44- We now see that the only parts of 
and that are needed for this calculation are those that correlate with (5^^^ or Thus all that 
is needed at the level of the fields is the basis So at the level of the fields and at two 

loop order we could write 

J(5) = + r(2)(5(2) + + T(3)<5® (B.5) 

to recover all the terms that enter in a two loop calculation. In the last line we reexpressed 5*^”^ in 
terms of . Besides the terms in the SPT two loop calculation the transfer functions will also pick 
up the counterterms to the extend that they are present in the data and correlate with first, second 
and third order fields. 
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